Optomechanical entanglement between a movable mirror and a cavity field 
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We show how stationary entanglement between an optical cavity field mode and a macroscopic 
vibrating mirror can be generated by means of radiation pressure. We also show how the generated 
optomechanical entanglement can be quantified and we suggest an experimental readout-scheme to 
fully characterize the entangled state. Surprisingly, such optomechanical entanglement is shown to 
persist for environment temperatures above 20K using state-of-the-art experimental parameters. 

PACS numbers: 03.67.Mn, 42.50.Lc, 05.40.Jc 



Entanglement, "t/ie characteristic trait of quantum 
mechanics" [l| , has raised widespread interest in different 
branches of physics. It provides insight into the funda- 
mental structure of physical reality [2] and it has become 
a basic resource for many quantum information process- 
ing schemes Q. So far entanglement has been exper- 
imentally prepared and manipulated using microscopic 
quantum systems such as photons, atoms and ions 
Nothing in the principles of quantum mechanics prevents 
macroscopic systems to attain entanglement. However, 
the answer to the question as to what extent entangle- 
ment should hold when going towards "classical" systems 
is yet unknown 5] . Therefore it is of crucial importance 
to investigate the possibilities to obtain entangled states 
of macroscopic systems [6] and to study the robustness 
of entanglement against temperature [7|. Experiments 
in this direction include single-particle interference of 
macro-molecules fsj, the demonstration of entanglement 
between collective spins of atomic ensembles Q, and of 
entanglement in Josephson-junction qubits [lol |. Mechan- 
ical oscillators are of particular interest since they re- 
semble a prototype of "classical" systems. Thanks to 
the fast-developing field of microfabrication, micro- or 
nano-mechanical oscillators can now be prepared and 
controlled to a very high precision llj. In addition. 



several theoretical proposals exist that suggest how to 



reach the quantum regime for such systems 12|. Ex- 
perimentally, quantum limited measurements have been 
developed that could allow ground state detection [isj . 
However, quantum effects in mechanical oscillators have 
not been demonstrated to date. 

Optomechanical coupling via radiation pressure d4] is 
a promising approach to prepare and manipulate quan- 
tum states of mechanical oscillators. Proposals range 
from the quantum state transfer from light to a me- 
chanical oscillator to entangling two such oscillators 
TsL liL 17, 18 1 . In this paper we propose an experimen- 



tal scheme to create and probe optomechanical entangle- 
ment between a light field and a mechanical oscillator. 
This is achieved using a bright laser field that resonates 
inside a cavity and couples to the position and momen- 
tum of a moving (micro)niirror. The proposal is based 
on feasible experimental parameters in accordance with 
current state of the art optics and microfabrication. In 
contrast to other proposals 15, [3] it neither requires 
non-classical states of light nor temperatures close to the 
oscillator's ground state. Entanglement is shown to per- 
sist above a temperature of 20K. We begin by modelling 
the system and its coupling to the environment by using 
the standard Langevin formalism. Then we solve the dy- 
namics and quantify the entanglement generated in the 
stationary state. Finally we discuss a suitable experimen- 
tal apparatus capable of measuring the entanglement. 

We consider an optical Fabry-Perot cavity in which 
one of the mirrors is much lighter than the other, so that 
it can move under the effect of the radiation pressure 
force. The motion of the mirror is described by the exci- 
tation of several degrees of freedom which have different 
resonant frequencies. However, a single frequency mode 
can be considered when a bandpass filter in the detection 
scheme is used [3l and mode-mode coupling is negligi- 
ble. Therefore we will consider a single mechanical mode 
of the mirror only, which can be modeled as an harmonic 
oscillator with frequency Wm- The Hamiltonian of the 
system reads [2^ 

Ti ~ hwco'a H T^{p^ + q^) — hGoo'aq (1) 
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where q and p {[q,p] — i) are the dimensionless posi- 
tion and momentum operators of the mirror, a and 



1) are the annihilation and creation operators 



of the cavity mode with frequency Wc and decay rate k, 
and Go = {wc/ L)yjh/mwm is the coupling coefficient. 
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with L the cavity length in the absence of the intracavity 
field and m the effective mass of the mechanical mode 
[l^. The last two terms in Eq. ([T]) describe the driving 
laser with frequency Wq and E is related to the input 
laser power P by \E\ = ^J^P^JhuiQ. 

A proper analysis of the system must include photon 
losses in the cavity and the Brownian noise acting on 
the mirror. This can be accomplished by considering the 
following set of nonlinear Langevin equations, written in 
the interaction picture with respect to hwoa'^a 



q = W,nP, 

p = -w,nq - -ymP + Goa^a + ^, 
a — —(K + iAQ)a + iGoaq + E 



(2a) 
(2b) 

2Ka"\ (2c) 



where Aq = Wc — wq and 7m is the mechanical damping 
rate. We have introduced the vacuum radiation input 
noise a*" , whose only nonzero correlation function is [21| 



(a™(Oa"^^(t')> ^6{t-t'), 



(3) 



and the Hcrmitian Brownian noise operator ^, with cor- 
relation function [231 
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[ks is the Boltzmann constant and T is the mirror tem- 
perature). We can always rewrite each Heisenberg oper- 
ator as a c-number steady state value plus an additional 
fluctuation operator with zero mean value, a = + (5a, 
q = qs + 5q, p = ps+ Sp. Inserting these expressions into 
the Langevin equations of Eqs. ([2]), these latter decouple 
into a set of nonlinear algebraic equations for the steady 
state values and a set of quantum Langevin equations for 
the fluctuation operators [23]. The steady state values 
are given hy Ps = 0, gs = Go|asp/u;,„, as = £'/(k + «A), 
where the latter equation is in fact a nonlinear equation 
determining the stationary intracavity field amplitude 
as, since the effective cavity detuning A, including radi- 
ation pressure effects, is given by A = Aq — Gl\as\'^ /wm- 
The parameter regime relevant for generating optome- 
chanical entanglement is that with a very large input 
power P, i.e., when \as\ 3> 1. In this case, one can safely 
neglect the nonlinear terms 6a^6a and 6a6q and gets the 
linearized Langevin equations 



Sq = WmSp: 

5p = -WmSq - ^mSp + GSX 

SX = -kSX + ASY + V2i^X"\ 
SY = -k6Y ~ ASX + G6q + V2^Y' 



(5a) 
(5b) 
(5c) 
(5d) 



where we have chosen the phase reference of the cav- 
ity field so that as is real, we have defined the cav- 
ity field quadratures 6X = {6a + 5a^)/y/2 and 5Y = 
{5a — 5a^)/i-\/2, and the corresponding Hcrmitian in- 



'^)/i^/2. What is relevant is that the quan- 



tum fluctuations of the fleld and the oscillator are now 
coupled by the much larger effective optomechanical cou- 
pling G = GqUsV^, so that the generation of significant 
optomechanical entanglement becomes possible. 

When the system is stable it reaches a unique steady 
state, independently of the initial condition. Since the 
quantum noises ^ and a™ are zero-mean quantum Gaus- 
sian noises and the dynamics is linearized, the quantum 
steady state for the fiuctuations is a zero-mean bipartite 
Gaussian state, fully characterized by its 4 x 4 corre- 
lation matrix Vij — {{ui{oo)uj{oo) + Uj{oo)ui{oo))) /2, 
where u^{oo) = (^g(oo), (5p(oo), (5X(oo), (5F(oo)) is the 
vector of continuous variables (CV) fluctuation opera- 
tors at the steady state {t oo). Defining the vector of 
noises n^{t) ^ {0,({t),V2^X'"{t),V2i^Y'"{t)) and the 
matrix 
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Eqs. ((5|) can be written in compact form as u{t) = Au{t) + 
n{t), whose solution is M(t) = M{t)u{0) + dsM{s)n{t — 
s), where M{t) = exp{Ai}. The system is stable and 
reaches its steady state when all the eigenvalues of A 
have negative real parts so that M(oo) = 0. The stability 
conditions can be derived by applying the Routh-Hurwitz 
criterion [2^ . yielding the following two nontrivial con- 
ditions on the system parameters. 



27„K [A^ + A2(7^ -f 27„K + 2k2 _ 2w: 
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WmG^A(7„ 
A2 -I- k2) - w„,G^A > 0, 



^) (7a) 
2k)2 > 0, 

(7b) 



put noise operators X* 



{a'' 



)/\/2 and Y'' 



which will be considered to be satisfied from now on. 
When the system is stable one gets 

/"OO />oo 

V^J^y, dsM,k{s)M,i{s')^ki{s-s'), {8) 

k,l -^0 

where $fei(s-s') = {{nk{s)ni{s') + ni{s')nk{s))) /2 is the 
matrix of the stationary noise correlation functions. Due 
to Eq. ([ll, the mirror Brownian noise ^(i) is not delta- 
correlated and therefore does not describe a Markovian 
process [22l|. However, quantum effects are achievable 
only using oscillators with a large mechanical quality fac- 
tor Q — Wmllm S> 1. In this limit, ^(i) becomes delta- 
correlated [251] ■ 

{mm + mm) /2 ^ i,n (2n + 1) 5{t - to, (9) 

where n = {explTiWrn/ksT} — 1) ^ is the mean ther- 
mal excitation number, and one recovers a Markovian 
process. As a consequence, and using the fact that 
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the three components of n{t) are uncorrelated, we get 
^kiis - s') = DkiS{s - s'), where D = Diag[0,7m(2n + 
1),K, k] is a diagonal matrix, and Eq. ^ becomes V = 
dsM{s)DM{sY . When the stabihty condit ions are 
satisfied, M(oo) = and one gets the following equation 
for the steady-state CM, 



AV + VA^ = -D. 



(10) 



Eq. (fTO|l is a linear equation for V and can be straight- 
forwardly solved, but the general exact expression is too 
cumbersome and will not be reported here. In order to 
establish the conditions under which the optical mode 
and the mirror vibrational mode are entangled we con- 
sider the logarithmic negativity Ej^, a quantity which 
has been already proposed as a measure of entanglement 
(2^ . In the CV case Ej^ can be defined as 27 1 



,1/2 



(11) 

1/2 



£'a^ — niax[0, — In 2?] ] , 

where 77- = 2-^/'^\t.{V) - [i:{Vf - 4,AeiV\ 

with S(V^) = dot A + det B — 2 det C, and we have used 
the 2x2 block form of the CM 



V 



A C 
B 



(12) 



Therefore, a Gaussian state is entangled if and only if 
r]^ < 1/2, which is equivalent to Simon's necessary 
and sufficient entanglement non-positive partial trans- 
pose criterion for Gaussian states [28| , which can be writ- 
ten as 4dety < E - 1/4. 

We have made a careful analysis in a wide parameter 
range and found a parameter region very close to that of 
recently performed optomechanical experiments (29| . for 
which a significative amount of entanglement is achiev- 
able. Fig. 1 shows Ej^ versus the normalized detuning 
A/wm, for two different masses, 5 and 50 ng: optome- 
chanical entanglement is present only within a finite in- 
terval of values of A around A ~ Wm- The robustness 
of such an entanglement with respect to the mirror's en- 
vironmental temperature is shown in Fig. 2. The rel- 
evant result is that for the 5 ng mirror optomechanical 
entanglement persists for temperatures above 20K, which 
is several orders of magnitude larger than the ground 
state temperature of the mechanical oscillator. For the 
50 ng mirror entanglement vanishes at lower tempera- 
tures (Fig. 2). Figs. 1-2 refer to Q = 10^ but we found 
that entanglement persists even for Q ~ lO'', although it 
becomes much less robust against temperature. In this 
case, entanglement persists up to 3(1) K for a 5(50) ng 
mirror. 

We finally discuss the experimental detection of the 
generated optomechanical entanglement. In order to 
measure Ej\f at the steady state, one has to measure all 
the ten independent entries of the correlation matrix V. 
This has been recently experimentally realized [sot for 




FIG. 1: Plot of the logarithmic negativity E// as a function 
of the normahzed detuning A/wm in the case of an optical 
cavity of length L = 1 mm, driven by a laser with wavelength 
810 nm and power P = 50 mW. The mechanical oscillator has 
a frequency Wm/Svr — 10 MHz, a damping rate yrn/2n = 100 
Hz, and its temperature is T = 400 mK. Full line refers to a 
mass m — 5 ng, and finesse J- = 1.07 x 10*, while the dashed 
line refers to a mass m = 50 ng and finesse = 3.4 x 10**. 
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FIG. 2: Plot of the logarithmic negativity Ej\f versus the 
mirror temperature. The full line refers to a mass m = 5 
ng, detuning A — Wm and finesse T = 1.07 x 10'*; the dashed 
line refers to a mass m = 50 ng, A — Wmfi, and finesse 
T = 3.4 X 10*. The other parameters are those of Fig. 1. 



the case of two entangled optical modes at the output of 
a parametric oscillator. In our case, the measurement of 
the field quadratures of the cavity mode can be straight- 
forwardly performed by homodyning the cavity output 
using a local oscillator with an appropriate phase. Mea- 
suring the mechanical mode is less straightforward. How- 
ever, if we consider a second Fabry-Perot cavity C2, ad- 
jacent to the first one and formed by the movable mirror 
and a third fixed mirror (see Fig. 3), it is possible to 
adjust the parameters of C2 so that both position and 
momentum of the mirror can be measured by homodyn- 
ing the C2 output. In fact, assuming that the movable 
mirror has unit reflectivity at both sides so that there is 
no light coupling the two cavities, the annihilation opera- 
tor of the second cavity, 02, obeys an equation analogous 
to the linearized version of Eq. 



6d2 



-{k2 + iA2)Sa2 + iG2a2Sq + ^2^4" W> (13) 



where K2, A2, a2, CL™(t) are the bandwidth, the ef- 
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fective detuning, the intracavity field amplitude, and 
the input noise of C2, respectively. Moreover, G2 = 
{w c2 / L2) y'fi / mwm, where Wc2 and L2 are the frequency 
and the length of C2. The presence of the second cavity 
affects the mirror dynamics, which is no more exactly de- 
scribed by Eqs. ([5]). However, if C2 is driven by a much 
weaker intracavity field so that \a2\ |as|, its back- 
action on the mechanical mode can be neglected and the 
relevant dynamics is still well described by Eqs. ([5]). 




Detector 



Loser 



FIG. 3: Schematic description of the proposed experiment, 
including the second Fabry-Perot cavity on the right for the 
detection of the mechanical motion. 

If we now choose parameters so that A2 = Wm 2> 
k2,G2\oi2\, we can rewrite Eq. (jl3p in the frame ro- 
tating at A2 — Wm for the slow variables 5d{t) = 
5o{t) exp{— iwmi} and neglect the terms fastly oscillat- 
ing at the frequency 2w,m so to get [l5l. [l8| 



G2^2 ~ 

5a2 — — K2<5a2 + i — -j^Sb 
\/2 



(14) 



where 5h — {iSp+Sq)/\/2. If K2 '3> G2|a2|/\/2, the cavity 
mode adiabatically follows the mirror dynamics and one 
has Sa2 — i{G2Ct2l K,2^)5b - 
A/2K2<5a2 — a™ ! finally get 



■ 727^02" 



{t). Using = 



it). 



(15) 



showing that, in the chosen parameter regime, the out- 
put light of C2 gives a direct measurement of the mirror 
dynamics. By changing the phases of the two local os- 
cillators and by measuring the correlations between the 
two cavity outputs one can determine all the entries of 
the CM V and from them numerically extract the loga- 
rithmic negativity Ej^ by means of Eq. pip . 

In conclusion, we have shown that a Fabry-Perot cavity 
with an oscillating micro-mirror and driven by coherent 
light can produce robust and stationary entanglement 
between the optical intracavity mode and the mechani- 
cal mode of the mirror. The amount of entanglement is 
quantified by the logarithmic negativity and surprisingly 
robust against increasing temperature: for experimental 
parameters close to those of recently performed experi- 
ments [2^ entanglement may persist above 20K in the 
case of a 5 ng mechanical oscillator. Finally, we suggest 
a readout scheme that allows a full experimental charac- 
terization of the CV Gaussian steady state of the system 
and hence a measurement of the generated entanglement. 
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